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Abstract 


Due to the exponential high gravitational red shift near the event horizon of a black hole, it 
might appear that the Hawking radiation would be highly sensitive to some unknown high energy 
physics. To study effects of any unknown physics at the Planck scale on the Hawking radiation, 
the dispersive field theory models have been proposed, which are variations of Unruh’s sonic black 
hole analogy. In this paper, we use the Hamilton-Jacobi method to investigate the dispersive field 
theory models. The preferred frame is the free-fall frame of the black hole. The dispersion relation 
adopted agrees with the relativistic one at low energy but is modified near the Planck mass m p . 
The corrections to the Hawking temperature are calculated for massive and charged particles to 
O (m~ 2 ) and neutral and massless particles with A = 0 to all orders. The Hawking temperature of 
radiation agrees with the standard one at the leading order. After the spectrum of radiation near 
the horizon is obtained, we use the brick wall model to compute the thermal entropy of a massless 
scalar field near the horizon of a 4D spherically symmetric black hole and a 2D one. Finally, the 
luminosity of a Schwarzschild black hole is calculated by using the geometric optics approximation. 
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I. INTRODUCTION 

Soon after Stephen Hawking demonstrated that quantum effects could allow black holes 
to radiate a thermal flux of quantum particles^, it was realized that there was the trans- 
Planckian problem with the calculation j^j]. Hawking radiation appears to come from the 
modes with huge initial frequencies, well beyond the Planck mass m p , which experience 
exponential high gravitational red-shifting near the horizon. So the Hawking radiation 
relies on the validity of quantum field theory in curved spacetime to arbitrary high energies. 
On the other hand, quantum field theory is considered more like an effective field theory of 
an underlying theory whose nature remains unknown. This observation poses the question 
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of whether any unknown physics at the Planck scale could strongly influence the Hawking 
radiation. 

To study the trans-Planckian problem, a hydrodynamic analogue of a black hole radiation 


was considered [3fl. Following the Unruh’s work, there have been a lot of studies to understand 

nn 

the dispersive held theory models|3H13|. which focused on studying the effect on the Hawking 
radiation due to modifications of the dispersion relations of matter fields at high energies. 
Similar to the original method for deriving the Hawking radiation, the energy fluxes for 
outgoing radiation were usually obtained by calculating the Bogoliubov transformations 
between the initial and final states of incoming and outgoing radiation. In most works, the 
Hawking effect could be recovered at leading order under some suitable assumptions, which 
have been briefly reviewed in 0Q. 

After the Hawking’s original derivation, there have been some other methods proposed to 
understand the Hawking radiation. Recently, a semiclassical method of modeling Hawking 


and attracted a lot of attention. This 
171 ], which is known as the null geodesic 


radiation as a tunneling process has been developec 
method was first proposed by Kraus and Wilczek|l6i. 
method. They employed the dynamical geometry approach to calculate the imaginary part 
of the action for the tunneling process of s-wave emission across the horizon and related it 
to the Hawking temperature. Later, the tunneling behaviors of particles were investigated 


using the Hamilton-Jacobi method 


18 


20]. In the Hamilton-Jacobi method, one ignores 


the self-gravitation of emitted particles and assumes that its action satisfies the relativistic 
Hamilton-Jacobi equation. The tunneling probability for the classically forbidden trajectory 
from inside to outside the horizon is obtained by using the Hamilton-Jacobi equation to 
calculate the imaginary part of the action for the tunneling process. Using the null geodesic 


method and Hamilton-Jacobi method, much fruit has been achieved 


21 


32]. Furthermore, 


the effects of quantum gravity on the Hawking radiation have been discussed in the Hamilton- 
Jacobi method. In fact, the minimal length deformed Hamilton-Jacobi equation for fermions 
in curved s pac etime have been introduced and the modified Hawking temperatures have been 


derived 


33 


38]. These have motivated us to use the Hamilton-Jacobi method to study the 


dispersive held theory models 


39]. In this paper, we focus on the dispersive models with the 


’ree-fall preferred fame, whereas those with the static preferred fame have been studied in 


39]. 


The remainder of our paper is organized as follows. In section QT| the deformed Hamilton- 
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Jacobi equations are derived for the dispersive models with the free-fall preferred frame. We 
then solve the deformed Hamilton-Jacobi equations to obtain tunneling rates for massive 
and charged particles to O (m“ 2 ) and massless and neutral particles to all orders. The 
thermal entropy of a massless scalar held near the horizon is computed in section IIIII using 
the brick wall model. In section IIV1 we calculate the luminosity of a Schwarzschild black 
hole with the mass M m p . Section [V] is devoted to our conclusion. Throughout the paper 
we take Geometrized units c — G — 1, where the Planck constant h is square of the Planck 
mass nip. 

II. DEFORMED HAMILTON-JACOBI METHOD 

In this section, we first derive the deformed Hamilton-Jacobi equation incorporating the 
modified dispersion relation (MDR) assuming that the preferred reference frame is the free- 
fall frame. We then solve the deformed Hamilton-Jacobi equation for the imaginary part 
of I which gives the tunneling rate T across the event horizon. We consider two cases, a 
massive and charged particle to O (m“ 2 ) and a neutral and massless particle with A = 0 to 
all orders. 

A. Deformed Hamilton-Jacobi Equation 

To study the deformed Hamilton-Jacobi method incorporating the MDR for the Hawking 
radiation, one first needs to choose the form of MDR in flat spacetime (the local free-fall 
frame) and generalizes it to curved spacetime. To be as general as possible, we will work 
with the MDR for a particle with mass m 

E 2 = F 2 ( p ) + m 2 , (1) 

where we define 

F(p)=pZ C ^’ ( 2 ) 

n=0 P 

nip is the Planck mass, Co = 1, and E and p are the energy and the norm of the momen¬ 
tum measured in some preferred reference frame, respectively. Note that the MDR (JTJ) is 
rotational invariant in 4D spacetime. To generalize the MDR ([I]) to curved spacetime with 
the metric we denote by tG the unit vector held tangent to the observers’ world lines, 
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which picks up a preferred frame. For a particle with the energy-momentum vector the 
energy E and the norm of the momentum p of the particle measured by these observers are 


E = p^, 
p 2 = E 2 - Pfjjf . 


(3) 

(4) 


The curved spacetime generalization of the MDR (|T|) with a preferred frame described by 
is obtained by plugging eqns. (J3J) and (J4]) into eqn. (JT|) . To obtain the deformed Hamilton- 
Jacobi equation incorporating the MDR, one needs to relate the classical action / to p^. 
In fact, it can be shown that, if / is a solution of the Hamilton-Jacobi equation, then the 
transformation equations give 

Pn = -dftl, (5) 

where — appears since p M = ( E , —p) in our metric signature. Replacing with / via 
eqn. (J5J) and putting eqns. (j3J) and (0J into eqn. (JT]) gives the deformed Hamilton-Jacobi 
equation. In the appendix of 39], the deformed Hamilton-Jacobi equation is also derived 
in the language of the effective field theory for a scalar field and a fermion one. There we 
considered a scalar/fermion with the mass m and charge q in a static black hole with the 
presence of electromagnetic potential A^. Neglecting self-interacting effective operators, we 
constructed the U (1) gauge invariant effective field theory incorporating the MDR. The de¬ 
formed Klcin-Gordon/Dirac equation was derived. The deformed Hamilton-Jacobi equation 
for scalars/fermions was then obtained using the WKB approximation. It was found there 
that the deformed scalar/fermionic Hamilton-Jacobi equation with respect to the preferred 
frame u fJ in the black hole background spacetime can be written as 


where 


T 2 = F 2 (X) + m 2 , 


T = -u> l («V + qA ,), X 2 = T 2 - («V + qA,f , 


( 6 ) 


(7) 


A^ is the black hole’s electromagnetic potential and q is the particle’s charge. 


As in 39], we here consider the black hole whose metric in the Schwarzschild-like coor¬ 
dinate is given by 


ds 2 = / (r) dt 2 — 


dr 2 

fir) 


— C (r 2 ) h a b (x) dx a dx b , 


( 8 ) 
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where / (r) has a simple zero at r = with /' (r/,) being finite and nonzero. The vanishing 
of / (r) at point r = r/, indicates the presence of an event horizon. We also assumed that 
the black hole is asymptotically flat which gives / (r) —» 1 as r —> oo. However, a more 
suitable coordinate for describing a specific family of freely falling observers is the Painleve- 
Gullstrand (PG) coordinate |5|, p, [40]. The PG coordinate anchored to the freely falling 
observers along the radial direction takes the form of 

ds 2 = dtp — [dr — v (r) dt p ] 2 — C (r 2 ) h ab ( x ) dx a dx b , (9) 

where v (r) is the velocity of the free fall observer with respect to the rest observer and t p 

measures proper time along them. The spacetime also has the event horizon at satisfying 
v (rh) = —1. We assume v < 0, dv/dr > 0 and v —> vq < 0 as r —> oo. Note that 

v < 0 means the infalling observers. For simplicity we specialize to the particular family 

of observers with no = 0 who start at infinity with a zero initial velocity. Since the vector 
field u M of the freely falling observers is tangent to the infalling world lines, one has for the 
infalling observers along the radial direction with v 0 = 0 that 


vf 1 = ^1, v (r), 0^ (10) 

in the PG coordinate and 

u 11 = ((ll) 

in the Schwarzschild-like coordinate. The fact that t p is the proper time along the infalling 
world lines means that u is equal to the gradient of t p , 


un = d^tp. ( 12 ) 

Using eqns. (HIT) and (fT2l) gives 

t p = t + J dr. (13) 

Substituting eqn. (TT3l) into the PG metric (J9]) and comparing it to the Schwarzschild-like 
coordinate (]8]), one finds 

v (r) = -y/ 1 - f(r). (14) 

In the PG coordinate, one can use eqn. (HU to show that eqn. (J7J) becomes 


T 
X 2 


- ( d t I + qA t ) - v (r) ( d r I + qA r ), 


( d r I + qA r ) + 


h ab (x) (d a I + qA a ) (d b I + qA b ) 
C(r 2 ) 


(15) 
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Note that dt and dt p are Killing vectors in the Schwarzschild-like coordinate and the PG 
coordinate, respectively. For the energy-momentum vector p eqn. (fT3j) gives us that its 
Killing energies associated with d t and d tp are the same. Let oo denote the Killing energies 
in the PG and Schwarzschild-like coordinates. Explicitly, one has u j = which 

is a constant. 

B. Massive and Charged Particle to O (jw^ 2 ) 

Since uo = —d tp I is the conserved energy of the particle, we can separate t from other 
variables. Thus, we employ the following ansatz for the action / 

/ = -ut + W (r) + 0 (x). (16) 

The vector potential A M is assumed to be given by 

A u = A (r) 5/rf, (17) 

which is true for charged static black holes in most cases. Putting the ansatz (1TB]) into eqn. 
dig), we have 


T — u r 
X 2 = p 2 r + 


r) — v (rj p r , 

h ab (x) d a O (x) dbO (x) 
C (r 2 ) 


(18) 


where Co (r) = co — qA t (r) and p r = d r W. The method of separation of variables gives the 
differential equation for 0 (x) 


h ab (x) d a Q (x) d b Q (x) = A, 


(19) 


where is A is a constant and determined by h ab (x). Thus, one has 

and eqn. (JBJ) becomes an ordinary differential equation for W (r). 

We consider a particle with mass m and charge q. Solving eqn. ([6]) for p r gives 

4 


( 20 ) 


± ±VA + \/l ~ f ~ ( \ ^ C i 

Pr = -7- W ( r ) =F —j 

/ rn- 


1 ± V(! - /) A w 3 (r) 


f 4 V A 


+ O (m; 4 ) , (21) 
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1 ' C(r'^) 

where +/— denotes the outgoing/ingoing solutions and A = 1-• Here, p+ has a 

pole at r = To obtain the residue of p+ at r = r^, one expands / (r) and C (r 2 ) at r = rh 

f (r) ~ 2k (r - r ft ) [l + r]K (r - r h ) + 0k 2 (r - r h f + pft 3 (r - r h ) 3 ] , 

C (r 2 ) ~ C (r^) [1 + ci/c (r - r h ) + c 2 k 2 (r - r h ) 2 + c 3 k 3 (r - r h ) 3 ] , 
a; (r) ~ d) (r h ) [1 + u\K (r - r h ) + a; 2 /t (r - r h ) + a; 3 K (r - r h ) \ 


( 22 ) 


where k — f () /2. Using the residue theory for the semi circles, we get 

CO (r h ) 7T 


Im W + (r) = 


K 


(1 + A), 


(23) 


where we define 
A — C ' 


SxX 


+ 5 m m 2 + <La; 2 (r/j) 


+ ^ K 4 ), 


2m 2 L C(rl) 

5\ = 3 + 2c 2 — 2c 2 + 12p + 12 ? 7 2 — 60 + c 3 (6 + 6p — 2aq) — 6uq — 6770 ;! + 2a; 2 , 

= — 1 + 12p 2 — 60 + p (4 — 6uq) — 2aq + 2 u 2 , 

S u = 40p 3 - 120 + 8p + p 2 (24 - 60^) - 3cc, + 2A6u l + 12a; 2 

— 2a; 2 + 12a; 2 — 12aqa; 2 + 77 (2 — 400 — 36aq + 24a; 2 + 24a; 2 ) — 6a; 3 . (24) 


On the other hand, one can use eqns. (j22j) to expand p~ at r = r h . It turns out that the 
residue of p~ at r = is zero. Hence, we have Im ITT (r) = 0. As shown in 39], the 
probability of a particle tunneling from inside to outside the horizon is 

2 


Pemit oc exp 


h 


(Im W+ - Im W-) 


(25) 


There is a Boltzmann factor in P em it with an effective temperature, which is 

T 0 


T eff = 


1 +A' 


(26) 


It is interesting to note that we have calculated A in the static preferred frame in 39]. For 
emitted particles with mass m and charge q, we found 


A = 


C 1 ! 


2m 2 L 


2 (3aq - p) cb 2 (r h ) + 


A 


C(r 2 h ) 


+ m 


+ o (m- 4 ) . 


(27) 
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C. Massless and Neutral Particle to All Orders 


We now work with a particle with m — 0 and q — 0. To get all order result, we let A = 0. 
Note that one has A = 0 for the solution in a 2D black hole or the s-wave solution in a 4D 
spherically symmetric black hole. In this case, we can use the following ansatz for the action 

/ 

/ = — cot + W (r ). (28) 

Hence, the deformed Hamilton-Jacobi equation becomes 


/ \2 2 
(U-Vp r ) =p r 


Y. 0 ' 


p 


2 n 


n , m 2n 
,n=0 H P 


(29) 


where p r = d r W. We will prove by induction that the solutions to eqn. take the form 
of 


pf = 


±o; 

1 ± v 


00 , ,2 i 

1 + 


c„ 




K U (i ± v) 


3 i-j 


where Ci j is determined by C n . In fact, it is easy to see that 


(30) 




p r = 


1 - 


cv 


i ± v m^(i ± i>y 
Suppose for some integer N > 1 such that 

±u 


, [3C? - c 2 (1 ± <01 , 0 


mt (1 ± v) 


Pr,N 


1 ± V 


N 2 i 

1 + 


a 


h3 


( 1±v ) 


3 i-j 


U 


m :: 


(31) 


(32) 


satisfy eqn. 


at O 


m* 1 
/ 2N-\-2 \ 

that the terms of O ( w 2 n +2 ) gives 
V m p J 


i . -I- rv^ 1 ( . . . 

• Plugging p r = p rN + ^ 2 iv +2 Txy into eqn. 

,p ±S-U 


± -Pjv-l (1 ± u) 
a = 


one finds 


(33) 


(1 ± n) 3(iV+1 ) ’ 

where Pat_i (x) is some polynomial of x with degree N — 1. This completes the proof that 


eqn. (T30l) are solutions to eqn. 
as 


to all orders. We define the residue of 


(i+W 


Res 


_ Rn 

(-1 j \ni r h J ) 

(1 +V) ) K 


at r = rh 


(34) 


where R\ = 1. Thus, one obtains 


LOTT 


Im W + = — (1 + A), Im W =0, 

K 


(35) 
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5iU 


2 i 


' i —1 


where we define 

A = S WUT With 5 i = ( Y1 C ^i+l-j 

i=l m P \j =0 

Using eqn. (|25|> gives the effective temperature 


T eff = 


1 + A 


(36) 


(37) 


D. Discussion 


When we use the residue theory for the semi circles to give eqns. 
sumption proposed in 


and fl35]h an as- 

is needed. The assumption requires that the singularity structure 


of d r I except the order of the pole at r = do not change after the MDR is introduced. It 
follows that oo < nip. A complete theory of quantum gravity might been needed to justify 
this assumption. 

In most works of the dispersive models, much attention have been paid to the modifica¬ 
tions of the asymptotic spectrum. Since the tunneling across the horizon takes place near 
the horizon, the near horizon spectrum of radiations is computed in the paper. The effects 
of scattering off the background need to be included if the asymptotic spectrum at infinity 
is considered. In most works, 2D spacetime has been considered. The higher order terms 
in the MDR violate conformal invariance of 2D spacetime, hence there is some scattering. 
Our calculations show that the spectrum of radiation near the horizon is close to a perfect 
thermal spectrum in the dispersive models. The thermal asymptotic spectrum has been 
recovered at leading order in previous studies. Thus, the energy fluxes of radiations are 
not significantly affected by the effects of scattering off the background. However, as noted 


m 


, the scattering effects might dramatically change the spectrum of radiations in the 


dispersive models with the static preferred frame. 


III. ENTROPY IN BRICK WALL MODEL 


Bekenstein and Hawking showed that the entropy of a black hole is proportional to the 
area of the horizon 


42 


44j . Although all the evidences suggest that the Bekenstein-Hawking 


entropy is the thermodynamic entropy, the statistical origin of the black hole entropy has 
not yet been fully understood. It appears that an unavoidable candidate for the statistical 


10 









origin is the entropy of the thermal atmosphere of the black hole. 

However, the entropy diverges when we attempt to calculate the entropy of the thermal 
atmosphere. There are two kinds of divergences. The first one is due to infinite volume of 
the system, which has to do with the contribution from the vacuum surrounding the system 
at large distances and is of little relevance here. The second one arises from the infinite 
volume of the deep throat region near the horizon. To regulate the divergences, t’ Hooft 45] 
proposed the brick wall model for a scalar field 0, where two brick wall cutoffs are introduced 
at some small distance r £ from the horizon and at a large distance L>ri, 


= 0 at r = rh . + r e and r = L. 


(38) 


In the following, we will use the brick wall model to calculate the entropy of a scalar held 
for a 4D spherically symmetric black hole and a 2D one. For a 4D spherically symmetric 
black hole, the entropy will be calculated to O (jn ~‘ 2 ). For the 2D black hole, we will obtain 
all order results in the cases with the static and free-fall preferred frames for comparison. 
For simplicity, we assume that the scalar held is massless and neutral. 

A. 4D Spherically Symmetric Black Hole 


For a 4D spherically symmetric black hole with the Schwarzschild-like coordinate 

dr 2 


ds 2 = / (r) dt 2 


f(r) 


C (r 2 ) ( dd 2 + sin 2 Odcf) 2 ) , 


(39) 


we have shown that A = (l + h 2 with the angular momentum l = 0,1, • • ■ and the 
corresponding degeneracy is 21 + 1 39]. Thus, the atmosphere entropy of a massless scalar 
held can be expressed in the form of 


s= f (2i+i) di f 


(40) 


where to is the Killing energy associated with t, l is the angular momentum, n (a;, l ) is the 
number of one-particle states not exceeding co with hxed value of angular momentum l , and 
s Ut i is the thermal entropy per mode. Taking the MDR corrections to both n (c u, l) and 
the Hawking temperature into consideration, we used the brick wall model to calculate this 
entropy to all orders in 39], where the static preferred frame is used to import a MDR to 
the black hole background. By contrast, here the MDR corrected atmosphere entropy of the 
black hole is computed to O 2 ) in the free-fall scenario. 
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For particles emitted in a wave mode labelled by energy co and l, we find that 


(Probability for a black hole to emit a particle in this mode) 
co 


= exp 


x (Probability for a black hole to absorb a particle in the same mode), 


T eff 

where T e ff is given by eqn. f[26j) . The above relation for the usual dispersion relation was 
obtained by Hartle and Hawking 50] using semiclassical analysis. Neglecting back-reaction, 
detailed balance condition requires that the ratio of the probability of having N particles 
in a particular mode to the probability of having N — 1 particles in the same mode is 


exp 


mode 


l eff 


. One then follows the argument in 


where we define 


n u ,i = n 


n {x) = 


co 


T eff 

1 


39] to get the average number n^j in the 

(41) 


(42) 


exp x — (—l) e 

Note that e = 0 for bosons and e = 1 for fermions. The von Neumann entropy for the mode 

is 

Su,i = K,j + (-1)1 In [1 + (-1 ) e n W}I ] - n w ,i In ra W)J . (43) 

Moreover, the entropy per mode s u .i can be put in the form of 

co 


Scj,1 — S 


eff 


(44) 


where the s (x) is given by 


s lx = 


-l) e exp x 


In 


exp x 


+ 


In [exp x — (—1)1 


expx —(—l) e |_ ex p re — (—1) 

Defining u = |r and expanding s Ut i to O (m~ : 2 ) give 

Slj,i ~ S (u) + s' ( u ) uA, 


exp a; 


-l) e 


(45) 


where one has 


A = 


C,S x z , CM 2 


+ 


2 m 2 p 


°u 2 . 


(46) 


(47) 


In the brick wall model, t’ Hooft’ found that the number of one-particle states not ex¬ 
ceeding uo fixed value l is 


n (c o, l ) = 


2ttH 


p r dr, 


(48) 
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where the integral ® p r dr was calculated in the Schwarzschild-like coordinate. Nevertheless, 


we calculate pf in the PG coordinate in section [TT] In 48, 491, the integral <j)p r dr has been 
found invariant under canonical transformations. Hence, the number states n (co, l) given in 
eqn. ()48l) is the same in the Schwarzschild-like and PG coordinates and one does not need 
to re-calculate it in different coordinates. Define the radial wave number k(r,l,co) by 


k± (r,l,u) = pf, 


± 


(49) 


as long as pf 2 > 0, and k ± (r, l, co) — 0 otherwise. With these two Dirichlet boundaries, one 
finds that the number of one-particle states not exceeding co fixed value l is 


n (c o, l ) = 


27 rh 


k + (r, l, co) dr + 


rrh+re 


k (■ r,l,uo)dr 


L Jr h +r e 


(50) 


The pf in eqn. 


held, one thus has 


are 


given by eqn. (l 2 T)l in section m For a massless and neutral scalar 


± ±VA + x/W C t 

Pr = -7- w =F — 

rrri 


1±v/A(1-/)1 4 ( 


co 


f 


H A 


+ O (rn p 4 ) , 


(51) 


where A = 1 — / 2 • 


Integrating by parts, one finds the entropy becomes 


S = - - / (21 + 1 ) dl / du 


ds U)[ 

du 


k (r, l, co) dr , 


(52) 


'r h +r e 


where we define k(r,l,co) = Plugging eqns. (15TT1 and (14511 into eqn. 

and performing the l integral which runs over the region where A > 0, we find that the 
entropy to O (r,n p 2 ) becomes 


5 


, C (rp r 0 3 
7 rh 3 
5GiTq u 2 s ( u ) 
2 mlf 3 


u 2 du 


r L C (r 2 ) dr 

Lr. WT) P 




2C,5 A C(r 2 )T 0 2 


m 


2 u 6 s' ( u) + ^ 9 ° w3g/ ( M ) 


3 fC (r 2 ) m2 


12 + 

5 + 2 32 


(53) 


where second and third terms in the bracket come from the MDR corrections to the Hawking 
temperature and the fourth term from the MDR corrections to n(co,l). Focusing on the 
divergent parts near horizon, we obtain for the nonnegative integers a and n 


f C a (r 2 ) dr 
L+r. C°(rl)f(r) 


n+1 


-i n— 1 rn,a 

k =0 


£n,a 

In K.r e , 

2 n+1 K 


(54) 
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oo 


where we expand / (r) n and ^ a |^ 2 | at 


r = r h 


f (r) n = 2 n K n (r - r h ) n ^/ '"kP (r - r h ) J , 


3=0 


C a (r 2 ) 

&Tr I) 


(r-r h )\ 


(55) 


t=0 


and define f£’ a = Ej= 0 f? +lc k-j ■ In eqn. 


, we neglect finite terms as nr £ —>■ 0 and 


terms involving L. Note that we define k = - y which is the surface gravity for the black 
hole and hence T 0 = Thus, the divergent part of entropy near the horizon to O (nip 2 ) is 

S ~ ^ ^ ' [ s (u) u 2 du f -(ci — 2 rf) In nr e 

l07T 4 J \KT e 

C (r 2 ) k 2 my f f 5C& ( /o 1 ’ 1 


167T 4 47T 2 

5CA 


s (u) vrdu 


KT r 


fi In kt £ 


6 

+ 5C\ 


1 f2,2 

J) 


dTTT ( Kr ^) fe 2 “ Z 2 ’ 2 ln 


2 - k 

y A 

32 ^ 4 - k 

k=0 


k =0 

3 ?4,1 




\ k —4 


O 2 X 3,1 f 1,1 

fc =0 


8 ray 


T/o’ 1 A 1 ’ 1 , 3/f 12/ 4 ’ 1 ' 


64 


+ 


32 


ln KTr 


(56) 


B. 2D Black Hole in Free-fall Scenario 


Consider a 2D black hole with the metric of 


ds 2 = f (r) dt 2 


dr 2 

/ M’ 


(57) 


in the Schwarzschild-likc coordinate. The atmosphere entropy of a massless scalar field is 

(58) 


f dn (w, () 
0=1 dOO -;- 


do; 

where is the entropy per mode. Define the radial wave number k (r, oo) by 

k ± (r,u)=pf, (59) 

as long as p^ 2 > 0, and k± (r, oo) — 0 otherwise. The pf in eqn. (l59|i are given in eqn. 
in section m The number of one-particle states not exceeding oo is 


n (oo) = 


2i:h 


rh+r e 


(k+ — k r ) dr. 


(60) 
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Defining the coefficients erf by 


2vr 


CD 


s=0 


2tt 


2s 


Aq= (mpKyi ( m P K Y S ..2S 


u 


one has for s u 


(61) 


s,, = s 


Ld 


T eff, 

00 r j(«) (n\ s,<l+ 2 r n q 


EE 

r =0 q=0 


s (q ’ ( u)u q+ - r a q r _ fm p K\ 2r 


q\ 


2vr 


where u = jr and we use eqn. (137)) for T e ff. Expanding (1 + v) 1 at r = r h 


(62) 


(1 + v )* = re* (r - r h y ^2 v )^ ( r _ r ^) J > 


j=o 


we hnd the divergent part of entropy near the horizon becomes 


S 


1 


oo oo r 


47T 2 


In Y. Y. 

r=0 i =0 q=0 


(2 i + 1) rjiG q (m p K\ 2i + 2r 


2n 


dus {q) (u) u q+2r+2i 


(63) 


_L y _j_ f nvn 2[,/;J1 yyy fEE ( ko) 2,+2r f dus M (u) 

4,r tT i(«r t )‘ V 2x ) iihii « ! ^ 2jt ) J 


where |~x~| is the smallest integer that is not less than x and we define 


2 — 1 


boo = 1 and r^o = 1, 

j=o 

3i+3[Z/3] + -Z 

= E ( 2s + 2 l'/3] + + 1) Ci 

J=0 


+[l/3\ + J V 3i-j-l+3[l/3] + ‘ 


(64) 


(65) 


C. 2D Black Hole in Static Scenario 


Following the conventions adopted in 39], the MDR for a massless scalar particle con¬ 
sidered here takes the form of 

P 2 = F 2 (E ), (66) 

where we define 


F(E) = Ey 


C„E 


2 n 


/yy~) 2n ^ 
71 = 0 rn P 


(67) 
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with Co = 1. For a particle with the energy-momentum vector p p , the energy uo and the 
norm of the momentum p of the particle measured by the static observers hovering above 
the 2D black hole with the metric (1571) are 

e = -£=, 

V7F) 

P 2 = f(r)p 2 r ■ (68) 

Relating p p to the action I by p p = —d p I gives the deformed Hamilton-Jacobi equation 

X 2 = F 2 (T), (69) 


where 


T = 


oo 


\/7M 


, X = \/7 {r)Pr, 


and oo = d^p p is the Killing energy with respect to t. Solving eqn. 
dehne the radial wave number k (r, oo) by 


(70) 


for p r , one could 


k{r,u) = \p r \, 


(71) 


as long as p 2 > 0, and k(r,oo) = 0 otherwise. The number of one-particle states not 
exceeding oo is 


" (w) = 


k r dr. 


The entropy per mode s u is 


Su, = s 


1'h~k r e 


00 


T f 


eff 


(72) 


(73) 


where s (x) is given in eqn. (l45|l and T e ff is the effective Hawking temperature. We calcu¬ 


lated T e ff in 


and it was given by 


T eff ~ 


1 +A’ 


where one has 


a = x>?<: 


,c o 


2k 


and p g fc and are defined in [39!]. Defining the coefficients a q s by 


2>k 07 . ) 

mi K 
k= l p 


(74) 


(75) 


A , = l'HVL\ 2, „2, 


2 7T 


u"' 1 <7? | —r— I u 


m p K\ 2s 2a 


s=0 


2vr 


(76) 
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we find the divergent part of entropy near the horizon becomes 


S 


1 v ' v ' v ' / 17lpK \ 2r, + 2r 2 Tl + 1 &r—q ^ 

—- c - 

n =0 r =0 < 7=0 


q\ 


en ,0 

J n—l 


hir, 



^ 1 («o)' 

where /^’ a with n> k > 0 are defined in eqn. (1541) . 


!1 2„ + 2r+, s (,) d% 


(77) 


D. Discussion 


In 


39| and this paper, we have calculated the divergent part of the near horizon at¬ 


mosphere entropy of a massless scalar held for a 4D spherically symmetric black hole in 
the static and free-fall scenarios, respectively. It appears that the divergent part in both 
scenarios can be presented in the form of a Laurent series with respect to r £ 


„0 °° 7^2 i 

~ — + Sq In nr e + ^ ° 


Si 


S ~ — + So In nr £ + > I > s) ( nr £ ) 3 + s l 0 In nr £ I , (78) 

*=i m P \j =i / 

where 5, = 2i + 1 in the static scenario and 5i = 3i + 1 in the free-fall scenario. Although 
we calculated the atmosphere entropy around the 4D black hole to O (m“ 2 ) in the free-fall 
scenario, the 2D black hole result suggests that eqn. (1781) might hold to all orders in this 
scenario. For s? and s{] in eqn. d78|) . we find that 

An 2 


«? = 


s o = 


720vr’ 

C (r l) k 2 

180 


(ci - 2 rj) , 


(79) 


where A = AirC (r^) is the horizon area. In the non-dispersive scenario ( m p — > oo), the terms 
and sjj In nr £ are the usual leading and subleading logarithmic divergent terms, respec¬ 
tively. Note that s® and have already been calculated in the non-dispersive scenario 45H47]. 

It seems from eqn. (l78l) that the near horizon divergence of the atmosphere entropy 
gets worse for the higher order corrections in the MDR as nr £ —> 0. However, the higher 
order contributions in eqn. ()78l) are always accompanied with the powers of the factor 
Tn 1 = ("IF") 2 ' Thus, one might hope that the higher order divergent problem would become 
less severe if r £ somehow can be related to rn p . One way to understand the value of r £ is 
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introducing the proper length for r £ as 


rrh+re 


£ = 


y/g^dr. 


(80) 


"Th 


The brick wall is put at r = + r £ to cut off the unknown quantum physics of gravity. In 

this sense, the invariant distance of the wall from the horizon e could be given by e ~ m p . 
Thus, we could define a such as e = am p . Indeed in the ’t Hooft’s original calculation, 


equating — to the Black hole’s Bekenstein-Hawking entropy Sbh — gives 

AvT* £ TTTlp 

nmi 


180tt 


(81) 


Note that eqn. (15U|) depends on the chosen coordinate system. In the scenario without the 
MDR, a natural choice is that £ is measured along a static time slice. Thus, eqn. (j80l) is 


calculated in the Schwarzschild-like coordinate 


45]. Assuming r £ <C r^, one finds from eqn. 


K£- 

~Y 


Thus, eqn. 


gives 


a = 


90tt’ 


(82) 


(83) 


where we reproduce’t Hooft’s result. 

In the static scenario, it is still natural to assume that e is measured along a static time 
slice. If we let £ = am p in eqn. (ITS]) , we find the atmosphere entropy around the horizon 
becomes 


5 


A 

4 m 2 v 


so + 2 sq In k m p + Finite terms as m p n —>■ 0, 


(84) 


where So was given in 39]. The leading divergent coefficient So is determined by the coef¬ 
ficients C n in the MDR (KJ71) and /j 1 and c“ which are dehned in eqn. (1551) . For a general 
black hole, /” and c“ could depend on the parameters of the black hole. However, they are 
pure numbers for a Schwarzschild black hole. Thus, for a Schwarzschild black hole, so does 
not depend on the black hole’s properties and the leading divergent term in eqn. (184|) scales 
with the horizon area A. 

In the free-fall scenario, one might prefer that the proper length £ is measured on a 
time slice orthogonal to the free-fall world lines 51]. In this case, eqn. (1801) for £ should 
be calculated in the PG coordinate and one then has £ = r £ . If one has £ = am p , the 
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atmosphere entropy around the horizon becomes 


S 


Ahi 


7207TOT7T, 


+ Sg In K'rrip + 


si 


(m p K) 


+ Finite terms as m p K —> 0, 


(85) 


where we define si = 


°l+2i 


jr. Moreover, eqn. (I55j) might suggest that effects of 


(27r) 2 W+ 2i 

i=max{l,Z—1} 

the MDR on the atmosphere entropy is nonperturbative in this case. Alternatively, inspired 
by the static scenario, one could choose r £ such that the higher order terms in eqn. (1751) 
have the same order of divergence as A-. Here, we could have 


r e = arripK 3 ’ 


where a is some constant. In this case, the atmosphere entropy becomes 

So 


S 


2 s° 


H-- In m p K + Finite terms as m p K —$■ 0, 


a (m p K) 3 


where we define 


s 0 — Sj 


i + E^ 


T ( 

i—1 1 


-2 i 


( 86 ) 

(87) 

( 88 ) 

For a Schwarzschild black hole, the terms in the square bracket in eqn. (188|) dose not 
depend on the black hole’s properties and the leading divergent term in eqn. (1571) scales 

4 _2 

with A/-c a nip 3 . In 51], the authors calculated the black hole horizon entanglement entropy 
for a massless scalar held with the MDR imposed in a free-fall frame. With the sub- or 
super-luminal dispersion with index n, they found that the entanglement entropy scales as 
An 1+ ™m v 


-i+ 


Following the argument proposed in 52], the authors in 53] obtained modified relations 
between the mass of a Schwarzschild black hole and its temperature and entropy. The argu¬ 
ment connecting a MDR and some modifications of the entropy of black holes is formulated 
in a scheme of analysis first introduced by Bekenstein 43]. In fact, for the MDR in eqn. (|T|) , 
the modified temperature of the black hole was given by 


T eff = —F 


4n \ 2 M 


m~ 


(89) 


where M is the mass of the black hole. The first law of black hole thermodynamics dSs = ^jr 
and eqn. (J2]) lead to the modified entropy of the black hole 

A 


S R = 


Artip 


+ 2ttCi In K'rrip + Finite terms as m p K —» 0, 


(90) 
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where A = 167 tM 2 and n = -Aj. It is interesting to note that the modifications of the 
entropy of black holes in eqn. (1901) are finite as m p n —y 0 except the logarithmic term. If 
one wants the same story for the atmosphere entropy obtained in eqn. fl78lh one conld have 
r £ = aK & ~ x nrip for some constant a where 0 < 5 < | in the static scenario and 0 < S < \ in 
the free-fall scenario. Hence, the atmosphere entropy becomes 

s° 

S ~-- —j + SQ^ln^mp + Finite terms as m p K —> 0, (91) 

a (ran p ) 

where the leading divergent term scales with AK 2 ~ s m~ s for a Schwarzschild black hole. The 
coefficient of the logarithmic term in eqn. (19U|) depends on C 1 , a coefficient of the MDR. 
However, we show that the coefficient of the subleading logarithmic term in the atmosphere 
entropy is irrelevant to coefficients of the MDR. It only depends on the position of the wall, 
r e and the properties of the black hole. 

For a 2D black hole with the Schwarzschild-like coordinate 

ds 2 = f (r) dt 2 - , (92) 

fir) 

the atmosphere entropy of a massless scalar can also be presented in the form of a Laurent 
series with respect to r £ 

00 \ 

s ~ In K'r E + + s 0 ln Kr ej , (93) 


where Si = 2 i in the static scenario and <5* = 3 i in the free-fall scenario. From eqns. (16411 
and (184(1 . one has Sq = — -A. In the static scenario, if we assume that the proper length £ 
is measured along a static time slice and £ = am p , the atmosphere entropy of a massless 
scalar becomes 


S ~ — ln Km n + Finite terms as Krn„ —> 0, 
6 p p 


(94) 


where the same leading logarithmic term was also obtained in 54] for the scenario without 


the MDR. In the free-fall scenario, if the proper length £ is assumed to be measured on a 
time slice orthogonal to the free-fall world lines and we let £ = am p , the atmosphere entropy 
of a massless scalar becomes 


1 1 00 

S ~-ln Km p + —- Y -— ^— r 

12 p 4 tH ^ rand’ 

n =1 v 


+ Finite terms as nm p —> 0, 


(95) 
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where ,s n can be determined by eqn. fIMjl . If we want that the modifications of the entropy 


of black holes in eqn. 


are finite, we could have r £ = an s l m d for some constant a where 


0 < 5 < | for some constant a. The entropy then becomes 


S ~ — — In Km p + Finite terms as nm p —>■ 0. 


(96) 


IV. BLACK HOLE EVAPORATION 


In the section, we discuss the MDR effects on the evaporations of a Schwarzschild black 
hole to O i 11 th e free-fall scenario. For simplicity, we assume that the emitted particles 
are massless. In 55], Page counted the number of modes per frequency interval with periodic 
boundary conditions in a large container around the black hole and divided it by the time 
it takes a particle to cross the container. He then related the expected number emitted per 
mode n U: i to the average emission rate per frequency interval -jp- by 


driu i duj 

= n u,i: 


(97) 


dt u ' , ‘27r/i’ 

for each mode and frequency interval (u,u + du). Following the same argument, we find 
that in the MDR case 

r] ' n ’ du dp r __ du 


UJ.l 


dt uj1 dpr 2nh u ’ l 27rh’ 

where is the radial velocity of the particle and the number of modes between the wavevec- 

OPr 

tor interval ( p r ,p r + dp r ) is jjr. Since each particle carries off the energy u, the total lumi¬ 
nosity is obtained from by multiplying by the energy u and summing up over all energy 


= n u ,r t 


u and l, 


i = E( 2 ' + D / 

i=o 




du 
2irh 


(99) 


However, some of the radiation emitted by the horizon might not be able to reach the 
asymptotic region. Before the radiation reaches the distant observer, they must pass the 
curved spacetime around the black hole horizon, which plays the role of a potential barrier. 
This effect on L can be described by a greybody factor from the scattering coefficients of 
the black hole. Actually, the greybody factor is given by |T) (co)| 2 , where Ti (u) represents 
the transmission coefficient of the black hole barrier which in general depends on the energy 
u and angular momentum l of the particle. Taking the greybody factor into account, we 
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find for the total luminosity 


L = (21 + 1) f \Ti(uj)\ 2 ujn U) i 

1=0 


du 
2irh 


( 100 ) 
e with the mass 


39] to compute 


The relevant radiation usually have the energy of order HM^ 1 for a black ho 
M, one hence needs to use the wave equations given in the appendix of 
|Xi(cu)| 2 accurately. However, solving the wave equations for |XJ(a;)| 2 could be very com¬ 
plicated. On the other hand, we can use the geometric optics approximation to estimate 
|Ti (w)| 2 . In the geometric optics approximation, we assume u M and high energy waves 
will be absorbed unless they are aimed away from the black hole. Hence |TJ (cu)| 2 = 1 for all 
the classically allowed energy oj and the angular momentum l, while |T, (E) | 2 = 0 otherwise. 
For the usual dispersion relation, the well-known Stefan’s law for black holes is obtained in 
this approximation. 

To find the classically allowed angular momentum l with fixed value of energy oj, we 
consider eqn. (JHJ) for a massless particle in the Schwarzschild black hole with the mass M. 
Solving eqn. (J6]) for A to O (m p 2 ) gives 


^ / \ 2 ^ 1 / \ 4 2 

= (CJ - Vp r ) - 5 - (CJ - Vp r ) - p r , 


, 2 


mt 


( 101 ) 


where we have A = (l + |) h 2 , v (r) = — \J^p- and C\ is given in eqn. (jTj). In the geometric 

optics approximation, p r is always a real number. In the non-dipsersive case ( m p —> oo), the 

^ 2 

maximum of the RHS of eqn. (11011) is U 2M . Thus, one has 


A < 


r 2 u 2 


( 102 ) 


l — 2M ’ 
r 

where the RHS has a minimum at r min = 3 M, which is 27 M 2 uj 2 . If the particles overcome the 
angular momentum barrier and get absorbed by the black hole, one must have A < 27 M 2 uj 2 . 
In the geometric optics_approximation, the Schwarzschild black hole is just like a black sphere 
of radius R = 3 3 / 2 Mjs3]. When the second term in the RHS of eqn. (11011) is included, the 


maximum of the RHS is shifted to 

H = 


or 


2CV 


or 


l _ 2M 


mt 


2M \ 4 ’ 


The minimum of r 2 Vt is then shifted to 


A m „ = 27 M 2 u 2 1 


(i-“) 


54C 2 C(; s 


mi 


(103) 
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Therefore, the particle must have A < A max to get absorbed by the black hole. The eqn. 
(110011 then becomes 


L = 


ojduj 
27 rh 3 


n 


to (1 + A) 


dX, 


(105) 


where eqn. (14TT) for n u _i is used. For a massless particle in the Schwarzschild black hole, eqn. 

(BID gi yes 


A = 


Ci 


2 


20 ok 


3A A 

4M2 ) ■ 


(106) 


Expanding n j to O (m“ 2 ) and then integrating eqn. (I105|l . we hnd for the emission 

of n s species of massless scalars and rif species of massless spin-1/2 fermions that the total 
lnminosity is 

7 


L = 


9 m 2 

40960 M 2 


n. 


mz 


+ 7 n f ) ~ (°- 7 3n s + 1-4171/) jTL + o ( 


mz 


M 3 


(107) 


To make a comparison with the static scenario, the results in 39] are given below for the 
form of MDR (Op) . In the static scenario, the correction to the Hawking temperature is 


A = 


Ci 

2rn 2 


4cC + 


A 


the maximum of the angular momentum is 


A ma ,v = 27 MV ( 1 - 


and the total lnminosity is 


L = 


9 mz 


409607 tM 2 


7 


AM 2 


6Ci07 5 


mi 


(108) 


(109) 


rri: 


n s + \ + C\ (0.48n s + 0.92-n/) jdL + Q 


m P 

M 3 


( 110 ) 


Note that the sign in front of Ci in eqn. (I107H is different from that in eqn. (II 10ft . For the 
sub-luminal dispersion relation with Ci < 0, it means that the total luminosity increases 
due to the MDR effects in the free-fall scenario while it decreases in the static scenario. 
In the geometric optics approximation, the black hole can be described as a black sphere 
for absorbing particles. The total lnminosity are determined by the radius of the black 
sphere R and the temperature of the black hole T. Note that one has R = y and 
T e ff « To (1 — A). In the static scenarios, the MDR effects increase the radius of the black 
sphere while they decrease the temperature of the black hole. The competition between the 
increased radius and the decreased temperature determines whether the luminosity would 
increase or decrease. It appears from eqn. (11101) that the effects of decreased temperature 
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win the competition. In the free-fall scenario, the MDR effects also increase the radius of 
the black sphere. Due to the minus sign in front of -AW in eqn. (I106jh the temperature of 
the black hole increases for A > |M 2 u 2 and decreases for A < ^M 2 (n 2 , but more slowly 
than in the static scenario. As a result, the effects of increased radius win the competition 
and hence the luminosity increases. The opposite story happens to the super-luminal case 
with Ci > 0. 

V. CONCLUSION 

In this paper, we used the Hamilton-Jacobi method to calculate tunneling rates of radi¬ 
ations across the horizon and the effective Hawking temperatures in the dispersive models 
with the free-fall preferred frame. After the near horizon spectrum of radiations was ob¬ 
tained, the thermal entropy of radiations near the horizon and luminosity of the black hole 
were computed. Our main results are: 

• In section QJ] we first derived the deformed Hamilton-Jacobi equations in the dispersive 
models with the free-fall preferred frame. The deformed Hamilton-Jacobi equations 
were then solved for d r I and the imaginary part of / was obtained. The corrections 
to the Hawking temperature were calculated for massive and charged particles to 
O {m ~‘ 2 ) and neutral and massless particles with A = 0 to all orders, respectively. It 
was found that corrections were suppressed by rn p . 

• In section m we used the brick wall model to compute the thermal entropy of a 
massless scalar field near the horizon of a 4D spherically symmetric black hole and a 
2D one. For a 4D black hole, the entropy near the horizon have been calculated to 
O (m^ 2 ) and could be written in the form of eqn. (178]) . The entropy became divergent 
as the wall approached the horizon. Various choices of the proper distance between 
the wall and the horizon and the corresponding entropies have been discussed. For a 
2D black hole, entropies in the static and free-fall scenarios have been calculated to all 
orders. The leading divergent terms were logarithmic. Nevertheless, their coefficients 
depended on choices of the proper distance between the wall and the horizon. 

• In section IIV1 we calculated luminosities of a Schwarzschild black hole with the mass 
M m p . We used the geometric optics approximation to estimate the effects of 
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scattering off the background. Comparison between the static scenario and the free- 
fall one has been given there. 
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